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The SU(3) r x Si7(3)£ linear sigma model is used to study 
the chiral symmetry restoring phase transition of QCD at 
nonzero temperature. The line of second order phase transi- 
tions separating the first order and smooth crossover regions 
is located in the plane of the strange and nonstrange quark 
masses. It is found that if the U(1)a symmetry is explic- 
itly broken by the U(1)a anomaly then there is a smooth 
crossover to the chirally symmetric phase for physical values 
of the quark masses. However, if the U(1)a anomaly is ab- 
sent, the region of first order phase transitions is significantly 
enlarged and it is found that there is a phase transition for 
physical values of the quark masses provided that the a me- 
son mass is at least 600 MeV. In both cases, the region of first 
order phase transitions in the quark mass plane is enlarged as 
the mass of the a meson is increased. 



The ultimate goal of relativistic heavy ion experiments 
is to probe the phase diagram of Quantum Chromody- 
namics (QCD). General theoretical considerations indi- 
cate that at sufficiently high temperatures there should 
be a transition from ordinary hadronic matter to a chi- 
rally symmetric plasma of quarks and gluons [Q. The 
order parameter for this phase transition is the quark- 
antiquark condensate. Results from lattice gauge theory 
predict the temperature of this transition to be about 
150 MeV The order of the phase transition, however, 
seems to depend very much on the number of quark fla- 
vors and their masses ||. 

Classically, the matter part of the QCD Lagrangian 
with Nf flavors is invariant under the symmetry group 
SU(N f ) r x SU(N f )t x U{1) A - The axial U(1) A symme- 
try is broken to Z{Nj)a by a nonvanishing topological 
susceptibility j| and the SU{Nj) T x SU{Nf)i symmetry 
is spontaneously broken to the diagonal group of vector 
transformations, SU(Nf) r+ e — SU(Nf)v, by a nonvan- 
ishing expectation value for the quark-antiquark conden- 
sate. The SU(N f ) r x SU(N f ) e x U(1) A group is also ex- 
plicitly broken by the effects of nonzero quark masses. It 
was shown by Pisarski and Wilczck that for three or more 
massless flavors, the phase transition for the restoration 
of the SU(Nf) r x SU(Nf)i symmetry is first order, while 
for two massless flavors the phase transition is second or- 
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The U(1)a symmetry may also be restored, if only par- 
tially, since instanton effects are Debye screened at high 
temperatures |^||. There are now two possibilities: ei- 
ther the U(X)a symmetry is restored at a temperature 
much greater than the SU(Nt) r x SU(Nf)g symmetry 
or the two symmetries are restored at (approximately) 
the same temperature 0. Recent lattice gauge theory 
computations have demonstrated a rapid decrease in the 
topological susceptibility at T c M and random matrix 
models also indicate that the two symmetries are restored 
simulaneously [^) . Perhaps more dramatically, it was also 
shown that the topological susceptibility vanishes at T c 
in the large- N c limit JToj] . On the other hand, the fate 
of the U(1)a anomaly in nature is not completely clear 
since instanton liquid model calculations indicate that 
the topological susceptibility is essentially unchanged at 
T c | p!l| . Additionally, other lattice computations which 
measure the chiral susceptibility find that the U(1)a sym- 
metry restoration is at or below the 15% level 12 13). 
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1. The phase diagram on the {m U} d,m a ) plane as ob- 
tained from lattice computations. These results are a com- 
pilation of data from the JLQCD and the Columbia groups 
taken from Refs. [3] and [14]. The plot is from Ref. [14]. 

Unlike the idealized massless quark limit, there are no 
general theoretical arguments which require that a phase 
transition exists for massive quarks. Indeed, some lat- 
tice simulations indicate that for physical quark masses, 
no phase transition occurs h JI^Ji^] . The general con- 
sensus from lattice computations is that in the plane 
of light quark masses (see Fig. [l]) there is a first order 
region bounded by a line of second order transitions. 
Outside this region, there is no phase transition, but 
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rather a crossover characterized by a rapid but smooth 
and continuous decrease of the quark-antiquark conden- 
sate. Given the present difficulties with performing lat- 
tice computations with realistic quark masses and a large 
number of sites, it is useful to complement the present 
lattice results with effective models that capture some 
of the relevant dynamics of QCD. Some work for three 
flavors has been done in this direction |l6|-]l9[]. In these 
works, the SU(3) r x SU(3)e linear sigma model was used 
to study the order of the chiral symmetry restoring phase 
transition as a function of the current quark masses with 
the ratio of the up-down to strange quark masses held 
fixed. In Refs. J^,[l7],[il| , a loop-expansion is used to 
compute the effective potential and in Ref. |H a mean- 
field analysis of this model was presented. The effects 
of the restoration of the U{1)a symmetry on the spec- 
trum of hadronic observables in heavy ion collisions was 
addressed within the context of this model in Ref. |^0| . 

In this paper, I present results concerning the or- 
der of the chiral symmetry restoring phase transition 
as a function of the current quark masses using the 
SU(3) r x SU(3)e linear sigma model without fixing the 
ratio of the masses. In addition, the effects of the U(1)a 
anomaly on the order of the QCD phase transition are in- 
vestigated. Here, the Cornwall- Jackiw-Tomboulis (CJT) 
pl| formalism is used to derive gap equations for the 
condensates and the tadpole-resummed scalar and pseu- 
doscalar nonet meson masses at nonzero temperature. 
The derivation of and the solutions to this set of equa- 
tions in a variety of limits has been presented elsewhere 
p2j. The results agree qualitatively with earlier stud- 
ies on a lattice |^,[l4|,[l5| and with other studies using 
the SU(3) r x SU{3) e linear sigma model In the 

presence of an explicit U(1)a symmetry breaking term, 
I find that for physical values of the strange and non- 
strange current quark masses, there is no phase transi- 
tion but rather a smooth crossover. For smaller values of 
the masses, the phase transition is first order with a line 
of second order transitions separating the first order and 
the crossover regions. In the absence of an explicit U(1)a 
symmetry breaking term, the region of phase transitions 
is greatly enlarged. In particular, if the a meson mass 
is greater than 600 MeV, then the transition is driven 
to first order for physical values of the quark masses. In 
both cases, the region of first order phase transitions is 
enlarged as the mass of the a meson is increased. 

The most general renormalizable theory compatible 
with the flavor symmetries of QCD is the SU(3) r x 
SU(3)i linear sigma model. While this model cannot 
account for the full dynamics of QCD, on the line of sec- 
ond order phase transitions the only relevant dynamics 
are determined by the symmetries of the theory. So, in 
the vicinity of this line, the use of the SU (3) r x SU (3)g 
linear sigma model is appropriate. Its Lagrangian is given 
by 



£($) = Tr (d^d^ - m 2 - A x [Tr ($ f $)] 2 
- A 2 Tr ($ f $) 2 + c [Dot ($) + Dct ($*)] 
+ Tr [#($ + $*)] . (1) 

Here, $ is a U(3) matrix defined by $ = T a (a a + i-K a ). 
The T a = A a /2 are the generators of U{3) where A a are 
the Gell-Mann matrices with Ao = \/2/3I. The T a are 
normalized such that Tr(T a Tb) = S a b/2. 

The parameters of the Lagrangian are the bare mass 
m, a background matrix field H = hoTo + hsTg, a cubic 
coupling c and two quartic couplings, Ai and A 2 . The 
various patterns of symmetry breaking and the parame- 
terizations of the coupling constants for this Lagrangian 
were studied in p2] | and will only be briefly reviewed 
here. For H_ = 0, c = and m 2 > 0, the Lagrangian has 
a global SU(3) r x SU{3) e x U(l) A symmetry. The ef- 
fects of the U(1)a symmetry breaking by a nonvanishing 
topological susceptibility (i.e. the presence of instantons 
in the QCD vacuum) are included by setting c ^= which 
reduces the symmetry to SU(3) r x SU(3)g. For nonzero 
H, chiral symmetry is explicitly broken. 

I assume that there are nonzero vacuum expectation 
values for the <Tq and a§ fields which I denote by tro and 
a a ■ After shifting these fields by their expectation values 
and following |23|| , the Lagrangian can be rewritten as: 

£ = ^[d^ad^CTa + df.TTad^TTa - (m|) Qb CT a CT b 

— (mp)ab^a^b} + [Qabc ~ 7j F abed &d) O a<JbO c 

4 

— 3 (Gabc + - HabcdOd) T^a^b^c — 2 H a bcd ^a^b^c^d 

— - Fabcd {o a ObO c ad + ^a^b^c^d) ~ h a O a , (2) 

where 

c 3 

Gabc — -^[d a bc ~ — (Saodobc + Sbod a Qc + ^cO^afco) 
2 

9 

+ -^doooSaoSboScf)] , 

Fabcd = -J (Sabred + SadSbc + S ac Sbd) 
A 2 

~l~ "7-7" (dabndncd ~t~ d a dndnbc d a cndnbd) j 
o 

"H-abcd — ~7~ SabScd ~1~ "7-7" (d a bndncd ~t~ facnfnbd 
4 o 

~H fbcnfnad) ; 
(m 2 s ) a b = m 2 Sab - 6 Gabc &c + 4 Tabcd O c O d , 
(mp) ab = m 2 Sab + 6 Gabc Oc + 4 Habcd O c O d . 

Here the summation runs over the index n only and d a bc 
and fabc are the symmetric and antisymmetric structure 
constants, respectively, of U(3). 

The a a fields are members of the scalar (J w = + ) 
nonet and the 7r a fields are members of the pseudoscalar 
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(J w = — ) nonet. The 7^,2,3 are the pions, the ^4,5,6,7 
are the kaons and the ttq and the 7rs are admixtures of 
the 77 and the 77' with mixing angle 6p. The situation 
with the scalar nonet is not as clear and still somewhat 
controversial Q. The cr and the a% are admixtures of 
the a and the /q(1370) with mixing angle 9s- The 01,2,3 
are identified with the ao(980) and the (74,5,6,7 with the 
k meson. 

The explicit symmetry breaking terms can be deter- 
mined (see, for instance, Ref. p^|), to be 

ho = ^=(mlU + 2m%f K ) (4a) 
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{ml U-m 2 K f K ) . 



(4b) 



The gap equations (Schwinger-Dyson equations) de- 
rived from the C JT effective potential plj in the tadpole- 
resummed approximation, or Hartree approximation, are 
found to be 



(Sabik))- 1 



-k +m 5 ab - QGabcVc + ^abcdCTc <Jd 



+ AT a bcd f S cd {k) + AHabcd I Vcd{k) , 
Jk Jk 

(Vab(k))- 1 = ~k 2 + m 2 5 ab + 6 Q ahc a c + m a i 
+ AHabcd f S cd {k)+AjF abcd [ V cd (k) 

Jk Jk 
h a = m 2 a a - 3 Gabc&b Oc - 3 Gabc J S c b(k) 

f 4 

+ 3Qabc J ~P cb {k) + -T abcd a d a b a c 



+ 4:T abcd <j d / S cb (k) + AH bcad a d / V cb {k) , (5) 

Jk Jk 

where in the last equation, a = 0, 8 and S ab (k) (V ab (k)) 
are the Green's functions for the scalar (pseudoscalar) 
mesons. Sos(k) and Vos(k), however, are nonzero on ac- 
count of the mixing between the singlet and the octet 
states. All other non-diagonal entries are identically 
zero. As such, it is necessary to rotate these Green's 
functions into the mass eigenbasis since only physical 
fluctuations can contribute to the masses: 

t^5 a6 (fc)t^=5i(fc)^ (6a) 

Uf a V ab {k)Uj b = V t {k)5 %3 , (6b) 

where U? a = S ia for i, a ^ 0, 8 and where U? a is given by 
an 0(2) rotation by Op in the 0-8 block. The definition 
for U° a is similarly given with Op — * 0g. The thermal 
integral arising from tadpole diagrams is 



Si{k) = 



d 3 k 



1 



(2rr) 3 ek[(M|)< 



exp ■ 



£k[(Mj)« 
T 



- 1 



and similarly for the pseudoscalar tadpole integrals, 
J k V cd (k). Here, e k [(Mj) 4 ] = (k 2 + (Af|)0 1/2 is the rel- 
ativistic energy of the ith scalar quasiparticle with mo- 
mentum k. I have neglected the vacuum contribution 



arising from the loop integrals. Implementing a system- 
atic renormalization scheme is difficult but possible in 
this approximation (see |p5|). The results, however, are 
not significantly altered. 

Since in the Hartree approximation, the gap equations 
do not have an explicit momentum dependence, we can 
assume that (S ab (k))^ 1 = — k 2 + M| where Ms depends 
on temperature but not momentum, and similarly for 
{Vabik))^ 1 ■ Equations (||) and (||) are then fixed point 
equations and can be numerically solved simultaneously 
as a function of temperature for M w , Mr, M v , M v i, M a , 
AI K , Mf a , M ao , ao, as, Op and 9s ■ The numerical solu- 
tions for a variety of parameters are given in Ref. [f2~2f . 

The condensate and mass gap equations are solved 
with fixed m, c, Ai and A2, while varying the background 
fields, ho and h%. The determination of the coupling con- 
stants is detailed in Ref. |22]| . For c ^ 0, the four cou- 
plings in the Lagrangian are fitted to yield the physical 
tree-level masses of the pion, kaon, <r, 77 and 77', while 
for c = 0, the remaining three couplings are determined 
from the physical tree-level masses of the pion, kaon, a 
and 77. The background fields are proportional to the 
current quark masses: m up = mdown = a{ho + h$/y/2), 
Estrange = b(ho — V^hg,) . For simplicity, I assume tem- 
perature independent proportionality constants, a and 
b. Requiring that m„ = 138 MeV, m K = 496 MeV, 
mup = "idown = 10 MeV and m stra nge = 150 MeV, gives 
a = 4.64 x 10" 6 [MeV]' 2 and b = 2.27 x 10~ 6 [MeV]~ 2 . 
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FIG. 2. The ao condensate for various values of the kaon 
mass with a pion mass of 100 MeV. The kaon mass is 80, 200 
and 300 MeV for curves (a), (b) and (c), respectively. Curve 
(a) is a first order phase transition, curve (b) is close to second 
order and curve (c) lies in the crossover region. 

To determine the order of the phase transition, I ex- 
amined the continuity of the order parameters as a func- 
tion of temperature. For a first order transition, the 
condensates are multivalued functions of temperature 
in the vicinity of the phase transition. For a smooth 
crossover, the condensates are smooth singlevalued func- 
tions of temperature and always nonzero. This behavior 
is demonstrated in Fig. ^. Only the nonstrange conden- 
sate is shown since both condensates exhibit qualitatively 
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the same behavior. 

The numerical results are plotted in Fig. ||. For 
c ^ 0, these results agree with those of lattice groups 

[IIHH. For m ° = 1000 MeV ' the auth ° rs of Ref - 

[ 18| report that the ratio of the critical current up- 
down quark mass to the physical up-down quark mass 
for m str angc/?7iup : down = 32 to be ~ 0.01. The corre- 
sponding value found in the present work is ~ 0.20. The 
larger value found here is most likely due to the inclusion 
of thermal fluctuations from the scalar and pseudoscalar 
nonets. For m up = mdown = 0, the transition is first or- 
der from zero strange quark mass to some critical strange 
quark mass. For m a = 800 MeV, this critical strange 
quark mass is about 16 MeV, while for m a — 900 MeV, 
it is 260 MeV. 
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FIG. 3. The lines of second order phase transitions in the 
plane of the nonstrange and strange current quark masses for 
m CT = 600 MeV, m CT = 800 MeV and m CT = 900 MeV. The 
cases where the U(1)a symmetry is explicitly broken by the 
axial anomaly, c ^ 0, are shown on the left, and the cases 
where the (7(1) a symmetry is exact, c = 0, are shown on 
the right. The physical mass point (m up ,down = 10 MeV and 
fJstrango = 150 MeV) is indicated by the diamond. 

For c = 0, however, the results are dramatically dif- 
ferent. In particular, the line of second order tran- 
sitions does not seem to approach the strange quark 
mass axis. For m a = 600 MeV, the physical point, 
m up = mdown — 10 MeV and m stra ngo — 150 MeV, is 
just outside the first order region. For larger values of 
the a meson mass, the physical point is well within the 
first order region. The results also seem to indicate that 
for c = there is a first order phase transition for three 
flavors provided only that one of the flavors is sufficiently 
heavier than the other two flavors. The departure of the 
second order phase transition line from the strange quark 
mass axis was also predicted using arguments from large- 



N c chiral perturbation theory in Ref. P6J. 

At this point, it should be mentioned that the Hartree 
approximation sometimes predicts a first order transition 
when the transition is actually second order. For exam- 
ple, renormalization group arguments predict a second 
order phase transition for the massless limit of the 0(4) 
linear sigma model jlj , while the Hartree approximation 
predicts a first order transition (see, for instance, Ref. 
p5[ ) . This is not a problem in the low quark mass region 
since the transition is expected to be first order. The lo- 
cation of the second order line should not be significantly 
affected. 

Additionally, the cubic and quartic couplings are fixed 
and temperature independent. The running of the cou- 
plings with temperature should be at most logarith- 
mic, while the integrals arising from the tadpole dia- 
grams depend quadratically on the temperature. So, 
it is reasonable that the running of the couplings does 
not qualitatively alter these results. On the other hand, 
if the Coleman- Weinberg mechanism is strongly opera- 
tive, some portion of the crossover region may actually 
be driven to first order p7[ . 
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